Information entropic superconducting microcooler 
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We consider a design for a cyclic microrefrigerator using a superconducting flux qubit. Adiabatic 
modulation of the flux combined with thermalization can be used to transfer energy from a lower 
temperature normal metal thin film resistor to another one at higher temperature. The frequency 
selectivity of photonic heat conduction is achieved by including the hot resistor as part of a high 
frequency LC resonator and the cold one as part of a low-frequency oscillator while keeping both 
circuits in the underdamped regime. We discuss the performance of the device in an experimentally 
realistic setting. This device illustrates the complementarity of information and thermodynamic 
entropy as the erasure of the quantum bit directly relates to the cooling of the resistor. 

PACS numbers: 74.50.+r,85.80.Fi,03.67.-a 



For the purpose of quantum computing, the coher- 
ence properties of superconducting quantum bits (qubits) 
should be optimized by decoupling them from all noise 
sources as well as possible. However, many interesting 
experiments can be envisioned also when the decoupling 
is far from perfect. One such experiment closely related 
to coherence optimization is using a qubit as a spectrom- 
eter 0, 0, Q for the environmental noise by monitoring 
the effect of the environment on the quantum two-level 
system. Here we focus on the opposite phenomenon, i.e. 
the effect of a qubit on the environment. Recently a 
superconducting flux qubit 4, 5] with a quite small tun- 
neling energy from the point of view of quantum com- 
puting was cooled using sideband cooling and a third 
level [6| from about 400 mK down to 3 mK. Motivated 
by this experiment we consider the possibility of using 
a single quantum bit as a cyclic refrigerator for environ- 
mental degrees of freedom. The utilized heat conduction 
mechanism is photonic which was recently studied also 
in experiment [7|. Besides the possible practical uses, 
the device is interesting physically as it directly illus- 
trates the connection between information entropy and 
thermo dynamical entropy. For related superconducting 
high-frequency cooler concepts see eg. Refs. [alSl- 

Here we study a flux qubit coupled inductively to two 
different loops shown in Fig. [IJl In loop j (j — 1, 2) we 
have a resistor Rj in series with an inductance Lj and 
a capacitance Cj . These form two damped harmonic os- 
cillators. The resistors are in general at different tem- 
peratures Xi and T 2 . The coupling of the qubit to these 
two admittances Y\ and Y 2 is assumed to be sufficiently 
large to dominate the relaxation of the qubit. This as- 
sumption can be easily validated by e.g. increasing the 
mutual inductance. The flux qubit is an otherwise su- 
perconducting loop except for three or four Josephson 
junctions with suitably picked parameters. In particu- 
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FIG. 1: (color online) Principle of the flux-qubit cooler, 
(a) Layout of the circuit, (b) Energy band diagram, (c) 
Schematic of the cooling cycle in the qubit temperature- 
entropy plane. 

lar one of the junctions is made smaller than others to 
form a two-level system. When biased close to half of the 
flux quantum <E>o = h/2e, the qubit can be described (in 
persistent current basis) by the Hamiltonian 

H/h = ~ (Aa x + ea z ) (1) 

where a x and a z are Pauli matrices, He — 2/ p ($ — $0/2) 
is the flux-tunable energy bias and <& is the controllable 
flux threading the qubit loop. Away from $ = $o/2 the 
eigenstates have the persistent currents ±/ p circulating 
in the loop. The tunneling energy HA results in an an- 



2 



ticrossing at $ = $0/2 and there the energy eigenstates 
do not carry average current. The resonant angular fre- 
quency of the qubit is ui = \/ e 1 + A 2 . 

Consider the ideal cycle shown in Fig. [TJd-c where the 
bias of the flux qubit is swept slowly (slower than A/2tt) 
between two extreme values E\ and £2 corresponding to 
two different energy level separations Yno\ and hw2- Let 
us further assume that ojj ~ ^LCj and Qj 3> 1, where 
LOLCj = 1/ y/LjCj and Qj = ^jLj/Cj/Rj. This choice 
guarantees that the qubit mainly couples to resistor Ri 
(R2) at bias point 1 (2). The cooling cycle consists of 
steps O, P, Q and R. First in step O the qubit has the 
angular frequency 0J2 and is allowed to thermalize. Be- 
cause of the bandwidth limitations imposed by the reac- 
tive elements, the qubit tends to thermalize with resistor 
R2 to temperature T2. In the next step P the flux bias is 
adiabatically changed to point 1 such that the level popu- 
lations do not change but the energy eigenstates do. The 
sweep is assumed to be however faster than relaxation. 
In point 1 the angular frequency is reduced to w\. Be- 
cause the level populations and therefore the Boltzmann 
factors do not change the qubit must now be at lower 
temperature T2 given by T2 = T2LO1/U2 in order to com- 
pensate for the change of the qubit splitting. Note that 
the quantum mechanical adiabaticity implies also ther- 
modynamical adiabaticity: while the energy eigenbasis 
changes the level populations and thus also entropy do 
not change. In step Q the qubit is allowed to thermalize 
to temperature T\ which results in heating of the qubit 
and in cooling of resistor 1 if T2 < T\. At this point 
the ideally pure quantum state of the qubit gets erased 
and information stored is lost. The entropy of the qubit 
increases, but locally the entropy of resistor 1 decreases 
such that one can say that some information is "stored" 
in the resistor as it cools but naturally with some loss. 
Finally in step R the qubit is adiabatically shifted back 
to frequency LO2 which results in heating of the qubit to 
the effective temperature 7\ = 7^2 /^i which is assumed 
to be higher than T2. The excess energy is dumped to 
admittance 2 when the cycle starts again from the be- 
ginning. Note that due to the condition T2 < T\ resistor 
1 can never be cooled below T2Wi/^2- Since there is no 
isothermal stage in the above cycle the present device 
is not even in principle a Carnot cooler but rather an 
Otto-type device. [lot 

The density matrix of the qubit with the resonant an- 
gular frequency u> at temperature T ((3 — (k-gT)^ 1 ) is 
given by 
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A £ \ f/3hu> 
—a x + —a z tanh — — 
uj oj I V 2 



(2) 



Using this the cooling power and the efficiency of the ideal 
cycle in Fig. [TJ; can be easily calculated. It is given by 
the area of the shaded region in the entropy-temperature 
plane below points P and Q. In principle one could solve 
for the effective temperature of the qubit along the line 



between points P and Q as a function of entropy given 
by S = — fcBTr(plnp). Alternatively, we can simply note 
that the expectation value of the energy stored in the 
qubit in point P is Ep = Tr(p cq (/32W2/wi, E\)Hi) while 
after relaxation we have Eq — Tr(/? eq (/3i,£i)ii/i), where 
Hi = H(e\) is the Hamiltonian at point 1. We thus get 
for the ideal cooling power 
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where / is the pump frequency. The cooling power 
achieves the maximum value of fuv\f /2 when the ther- 
mal population in step O (and P) is small and when 
the population in step Q is large, i.e. when $2^2 3> 1 
and Pifavi <C 1. Naturally a practical device has to be 
designed to fulfill the first condition always, in which 
case the smallest achievable temperature is on the or- 
der of fkvi/kB below which the cooling power decreases 
rapidly. The dynamic range could be made wider by a 
tunable A which can be achieved by splitting the small- 
est junction into a dc SQUID geometry. Another figure 
of merit is the ratio 77 of the heat removed from resis- 
tor 1 divided by the heat added to resistor 2. It can 
be obtained as the ratio of the shaded area divided by 
the sum of the hatched area and the shaded area, i.e., 
77 = (Eq-E p )/(E r -Eo) where E = Tr(p cq (/3 2 , £2)^2) 
and Er — Tr(p cq (/3iu>i/uj2, £2)^2). This simplifies neatly 
to 77 = uji/lu2 < 1 which is in harmony with the second 
law of thermodynamics. 

For more quantitative analysis we have to consider the 
details of the relaxation rates due to the baths. The 
Golden Rule transition rates due to resistor j are given 
by 
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where the positive sign corresponds to relaxation. The 
total thermalization rate is r^ h = + Tj. Here the 
unsymmetrized noise spectrum is given by 
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where ReYj (w) = R-'/[l + Q 2 { 
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J -^J-f \ is the real 
part of admittance of circuit j. The total relaxation rate 
is thus 
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To model the behavior of the device we utilize the Bloch 
master equation (see e.g Ref. (llj ) given in our case by 



M 



-BxM-TKM^-Mt,) 



-r t 2 h (M r 



■M T 



T 2 M ± , 
(7) 

where M = Ti(ap) is the "magnetization" of the qubit, 
and B = Ax + ez is the fictitious magnetic field. Note 
however that the z-component of B and M do correspond 
to real magnetic field and magnetization, respectively. In 
Eq. ([7]) Mil and M± are the components of the magnetiza- 
tion parallel and perpendicular to B, respectively. These 
are explicitly 



My = ^(AM X + eM z )(Ax + ez) 

A 2 M z ~AeM x 



M± = E2M 'Z? eM ' Z+M y y 



(8) 
(9) 
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Here Mt stands for the £-dcpcndcnt equilibrium mag- 
netization of a qubit at temperature T given explicitly 
by 



M T = 
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and T 2 = (rj h + r 2 h )/2 + r v is the dephasing rate. The 
possibility of pure dephasing at the rate T v has been 
included. In the simulation we neglect pure dephasing 
due to the intentionally large dominating thermalization 
rate. Equation ([7]) describes relaxation towards instan- 
taneous equilibrium with two competing rates due to 
two different thermal baths. Equations of this type are 
usually used in the stationary case, but for driving fre- 
quencies slower than A/K it should be also valid. As 
is obvious from Eq. 0, the qubit actually tends to re- 
lax towards an effective e-dependent equilibrium mag- 
netization (rl h M Tl +r t 2 h M T2 )/(r t 1 
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T 2 h ) at the rate 
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To illustrate the practical potential of the device we 
show in Fig. [2] the simulated cooling power with si- 
nusoidal driving of e(t) compared to the ideal case 
along with the actual loop in the entropy temperature 
plane. The heat flow Pj from resistor j to the qubit 
is simply obtained by integrating the product of the 
thermalization rate and the energy deficit, i.e., Pj = 
ffo /f dtT L [Tr(Pc q (&, e(t))H) - Tr(p(t)H)] . The den- 
sity matrix p(t) = ^M(t) ■ a is solved numerically using 
the Bloch equation (system is followed over a few periods 
until it has converged to the limit cycle) . We see that the 
actual simulated behavior does not significantly deviate 
at low / from the ideal behavior and that cooling pow- 
ers on the order of fW can be achieved with reasonable 
sample parameters. The oscillatory behavior at high / is 
interpreted as Landau-Zener interference 12, 13| . 

However, the cooling power has to be compared with 
realistic heat loads to evaluate the utility of the flux qubit 
cooler. On one hand, resistor 1 is subject to heat load 



FIG. 2: (color online) Example of the simulated cooling power 
with uji/2-k = A/2n = 5 GHz (e = GHz), u 2 /2tt = 20.62 
GHz (e = 20 GHz), Q 1 = Q 2 = 10, lo, = lo LCj and 
2(I p AM j ) 2 /(R j hcu j ) = 20 x 10 9 s~\ This can be achieved 
e.g. with J p = 200 nA, Mi = 29 pH, M 2 = 59 pH and 
R% = i?2 = 1 fi. The driving is sinusoidal, (a) The solid 
line illustrates the path in the T — S plane for the ideal cy- 
cle described in the text while the dashed (dotted) line is a 
result of simulation for / = 0.05 GHz (/ = 1 GHz) with 
Ti = T 2 = 0.3 x huj 2 /k B w 300 mK. (b) Simulated cooling 
power vs. / for the same temperatures as in (a) is shown 
with the dashed line while the solid line is the ideal result of 
Eq.|3](c-d) Same as (a-b) but with Ti = 0.5 x T 2 at 150 mK. 
The cooling threshold at 0.14 GHz in (d) is caused by finite 
Q-factor. 



from the phonons of the substrate on which the device 
rests. On the other hand, resistor 2 should be coupled 
well enough to phonon bath such that the unavoidable 
work done on it does not raise T 2 excessively. The heat 
flow between the electron system of resistor j and the 
phonon system is given by P e i- P h = — T^ h ) where 

Vj is the volume of resistor j and S is typically on the 
order of 10 9 Wm~ 3 K~ 5 . Thus resistor 1 needs to have a 
sufficiently small volume while resistor 2 should be large 
enough physically in order to serve as a heat sink. In ad- 
dition the photonic heat conduction between the resistors 
due to temperature gradient may in principle contribute 
also. Following an analysis similar to Ref. 14|, the heat 
flow from admittance ^2(0;) to Y\(lo) can be written as 



P-y = 



^ \4hto 3 M 2 ReY!(Lj)ReY 2 (Lu)(n 2 {uj) -ni(w))l 
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(11) 

where rij(u)) — [exp(f3jhuj — 1)] 1 are the boson occu- 
pation factors and M is the mutual inductance between 
the loops. For detuned high-Q resonators the photonic 
heat conduction turns out to be quite negligible. For 
instance for the values of Fig. [3] with M = 5 pH and 
fli = R 2 = 1 CI we get only P 1 = 2 x 10" 18 W even 



4 



0.3 



0.25 



0.2 



-0.15 



0.1 



0.05 



0.5 



f (GHz) 



1.5 



FIG. 3: (color online) Equilibrium temperature as a function 
of pump frequency for three different phonon bath temper- 
atures. The temperature of resistor 1 (volume 10~ 21 m 3 ) is 
shown with dashed line while the temperature of resistor 2 
(volume 10 -18 m 3 ) is shown with solid line. The bath tem- 
peratures T p h ~ T2 from top to bottom are 0.3 x fiu^/fcrj, 
0.2 x hu)2/k& and 0.1 x hu>2/kB- Otherwise the parameters 
are like in Fig. [21 



if T x = K and T 2 = 300 mK. Figure H illustrates the 
calculated equilibrium temperature versus operation fre- 
quency obtained numerically by finding the balance be- 
tween the dominating phononic heat conduction and the 
integrated cooling power. We see that almost a factor of 
2 reduction of T± is possible with realistic parameters. 

In practice the drop of T\ can be measured e.g. using 
an additional SINIS thermometer, in which resistor 1 will 
serve as the normal metal N. Its reading is sensitive to 
the electronic temperature of N only, and self-heating can 
be made very small. The resistors should be made out of 
thin film normal metal such as copper or gold with typ- 
ically sub 1 Q, square resistance. Volume can be picked 
freely. To get the resonant frequencies and quality factor 
as above we need L\ — 320 pH, C\ — 3.2 pF, L 2 = 80 
pH and C2 = 0.8 pF which are also realistic. For the 
inductor one may use either Josephson or the kinetic in- 
ductance of superconducting wire while the capacitance 
values are similar to those in typical flux qubits Q. To 
satisfy the conditions of the above numerical example we 
need quite large mutual inductances which however can 
be easily achieved using e.g. kinetic inductance [lij ]. The 
strong driving requires also rather large inductance be- 
tween the microwave line and the qubit, which should not 
result in uncontrolled relaxation. For instance, M mw =5 
pH coupling to the control line is acceptable as it would 
result in at most 3 x 10 7 s _1 relaxation rate assuming a 



50 environment at 0.3 K. This choice will not degrade 
the performance of the device significantly since driving 
is much faster. Yet sufficiently strong driving can be 
achieved with a modest 3 ac current. Fabrication 
process will require most likely three lithography steps. 

In conclusion, we have described a method of using a 
superconducting flux qubit driven strongly at microwave 
frequency to cool an external metal resistor. Here we con- 
sidered LC resonators to achieve the required frequency 
selectivity but a coplanar wave-guide resonator or a me- 
chanical oscillator could be used in principle, too. We 
demonstrated by a numerical example that it is possible 
to observe the associated temperature decrease experi- 
mentally. This effect is directly related to the loss of 
information and thus to the increase of entropy of the 
quantum bit. 
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